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We study the effect of extra dimensions on the process of massive Dirac fermion emission in the 
spacetime of (4 + n)-dimensional black hole, by examining the Dirac operator in arbitrary spacetime 
dimension. We comment on both bulk and brane emission and find absorption cross section and 
. . , luminosity of Hawking radiation in the low-energy approximation. 
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^-^ • I. INTRODUCTION 

' Scientists have devoted many years to the ongoing quest to unify the forces in Nature. Higher-dimensional theories 
I provide a promising framework for the unification of gravitation with other fundamental forces. In this context 
■ braneworld models Q with large extra dimensions point us the way out of the long standing hierarchy problem by 
. lowering the fundamental scale of gravity down to order of TeV. It has been also argued that mini black holes might 
(— I ' be created through high-energy particle collisions at TeV-energy scales. These objects are expected to evaporate 
•'^ , through Hawking radiation both in the bulk as well as on the brane. Mini black holes created in high-energy collisions 
^ • will undergo a number of phases, i.e. balding phase when black hole will emit mainly gravitational radiation, spin- 
O ] down phase during which black hole will loose its angular momentum through emission of Hawking radiation and 
^ Schwarzschild phase in which black hole will lose its actual mass by Hawking radiation. Finally in Planck phase 
quantum gravity theory is needed to study its behaviour. 
y—( , The TeV scale gravity opens up the possibility of producing black holes and observing their decay product. The 
I> ' aforementioned range of energy will be soon achieved by CERN Large Hadron Collider. One hopes that it proves or 
, merely restricts the parameter range (e.g., number or size of extra dimensions) of higher dimensional theories. 

Studies of particle emission from multidimensional black hole have their own long history. Namely, in Ref.Q 
massless scalar emission were studied in the spacetime of (4 -I- n)-dimensional Schwarzschild black hole while the case 
of massless spinor and gauge particles was treated in Q . Then, radiation emitted from higher-dimensional black holes 
] were considered both analytically and numerically (see, e.g., Refs.0-[^ for a non-exhaustive sampling of this widely 
treated subject). On the other hand, graviton emission in the bulk from a higher dimensional Schwarzschild black 
hole was elaborated in Refs. 0] , where it was established that the low-energy emission rate decreases with the number 
of extra-dimensions as was previously found for the case of bulk massless scalar field. 

The complexit y of the aforementioned problem in the background of a rotating (4 + n)-dimensional black hole was 
revealed in Refs.[lO|. The argument of the effect of mass on the emission spectrum in four-dimensions was quoted 
in [ll[ (see also where the particle and light motion in the vicinity of five-dimensional rotating black hole was 
' investigated). In Ref.[l^ the scalar massless Hawking emission into the bulk by a higher dimensional rotating black 
hole with only one non-zero angular momentum component in a plane parallel to the brane was investigated. It was 
found that the bulk emission remained significantly smaller than the brane one. It turned out also that the angular 
momentum loss rate had smaller value in the bulk than on a brane. Fermion emission in (4 -|- n)-dimensional rotating 
background obtained by projecting the higher dimensional black hole line element on the brane was studied in |14| 
where was revealed that any increase of dimensionality or angular momentum parameter of the black hole significantly 
enhanced all emission rates. The same conclusions were drawn in Ref.[l^. It was shown in [16;], that for spacetime 
of dimensionality greater than five bulk fermion emission dominated brane localized one, contrary to the conjecture 
that black holes radiated mainly on the brane. Recently, greybody factor of Schwarzschild black hole localized on a 
three-brane of finite tension, both in the low and high frequencies regimes were studied [l7j . 

The absorption probabilities of massless scalar field were also studied in more exotic backgrounds, e.g., in the 
spacetime of a rotating Godel black hole in minimal five-dimensional gauged supergravity as well as in the space- 
time of rotating Kaluza-Klein black hole with squashed horizon [l^. The Hawking radiation of (4 + n)-dimensional 
Schwarzschild black hole imbedded in de Sitter spacetime was investigated in Ref. [20| . The problem of massless scalar 
emission in the background of n-dimensional static black hole surrounded by quintessence was elaborated in Ref. plj . 
On the other hand, the emission spectrum of minimally coupled particles with spin which were Hawking radiated 
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from four-dimensional black hole in string theory was studied in [22l | , while the low-energy cross section for minimally 
coupled massless fermions was provided in Ref.[23|. The low-energy absorption cross section for massive fermions in 
the Schwarzschild background defined on the bulk was investigated in [2J| using Dirac equation in the traditional 
form. 

Our present work will be devoted to (4 + n)-dimensional black holes evaporating massive Dirac degrees of freedom. 
In our research we use the much easier method of solving massive Dirac equation than the traditional one. We 
shall restrict our attention to the static spherically symmetric case. It turned out that the treatment of fermions in 
spherically symmetric background may be simplified by using a few basic properties of the Dirac operator. These 
features enable us to find the second order differential equation which will be crucial in further investigations. In Sec. II 
we recall the basic features of Dirac operator which enable us to treat the aforementioned problem of massive Dirac 
fermions in arbitrary spacetime dimension. In Sec. Ill we shall focus on the low-energy regime and solve analytically 
Dirac massive fermion field equations by means of the matching technique combining the far field and near event horizon 
solutions. We find the analytical expression for the absorption probability and luminosity of Hawking radiation. We 
shall comment both on bulk and brane fermion emissions. Our analytical considerations are supplemented by plots 
expressing the dependences of the found quantities on various parameters of the considered spacetimes. Finally, in 
Sec. IV we shall state our conclusions. 



II. PROPERTIES OF THE MASSIVE DIRAC EQUATION 

In this section we recall the properties of Dirac operator (for the readers convenience) which enable us to simplify 
studies of the Hawking massive fermion radiation in higher dimensional black hole backgrounds. It happened that 
the treatment of fermions in spherical spacetime may be greatly simplified due to the few properties of the Dirac 
equation of motion. This attitude was widel y u sed in studies of the late-time behaviour of massive Dirac hair in the 
spacetimes of various black holes (see, e.g., |25l|-[27j). One should remark that our considerations are conducted in 
arbitrary spacetime dimension. 

The massive Dirac equation in a curved background may be written in the form as (for the convention we used see 
Refs-Hi,!!!) 



7''V^-mjV = 0, (1) 

where is the covariant derivative = 5^ + juj'^^^alb, M ^-^d a are tangent and spacetime indices. They are 
related by e° with a basis of orthonormal one-forms. The quantity uj"^^ = uj"-'' is associated with connection one-forms 
satisfying de° -I- A = 0, while the gamma matrices fulfill the relation {7", 7^} = 277"''. Now, we shall recall some 
basic properties of the Dirac operator^ = 7^V^ on an n-dimensional manifold. 

In what follows we assume that the metric of the underlying spacetime may be rewritten as a product of the form 

g^.^^dx^'dx" = gab{x)dx''dx'' + g„,n{y)dy"'dy"- . (2) 

The above metric decomposition will be subject to the direct sum of the Dirac operator, namely one obtains 

lp=lp^+lpy. (3) 

By virtue of a Weyl conformal rescaling defined by the following: 

9f,u = ^"^g^iu, (4) 

where is a conformal factor, the above considerations consequently provide that one gets 

p^|J = f2-^("+i)|)V^ , = n-i'^''-^^. (5) 

Having in mind that spherically symmetric form of the line element provides also conformal flatness for a static metric, 
one obtains 

ds^ = -A^dt^ + B^dr^ + C^dnl^2 , (6) 

where A — A(r), B = B(r), C = C(r) are functions only of the radial variable r, and the transverse metric rffJ^^j is 
independent on t and on r-coordinates. 
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Suppose then, that ^' is a spinor eigenfunction on the {n + 2)-diniensional transverse manifold fi. It leads to the 
relation 

= A*. (7) 
Using the properties given above one may also assume that the following is satisfied: 

P^P = mi/'. (8) 

It enables us to set the form of the spinor i/;, i.e., 

V' = -^^^X'^*- (9) 

Next, the explicit calculations reveal 

{l'dt+j'dy)x = Aim~^)x, (10) 

where we have introduced the radial optical distance (i.e., the Regge- Wheeler radial coordinate) dy — B/Adr and 
7°, 7^ satisfy the Clifford algebra in two spacetime dimensions. 

We remark that an identical result may be obtained if a Yang-Mills gauge field is present on the transverse manifold 
E. The only difference is that 



* = A*, (11) 



where is the Dirac operator twisted by the the connection A^. Thus, in the picture under consideration, 

assuming that ij) oc e~*"* one achieves the second order equation of the form as follows: 

III. GREYBODY FACTOR IN THE LOW-ENERGY REGIME 

In this section we shall concentrate our attention on finding the absorption probability in the low-energy regime. 
Greybody factors enable us to study the near horizon structure of black holes. This is of a great importance from 
the experimental point of view due to the fact that they modify the spectrum in the region of particle production. In 
general the spectrum of emitted particles depends on various factors such as spin of the particles, whether the particle 
is localized on brane or can propagate in the bulk. The greybody factor can be computed by finding the absorption 
cross section for the type of particle, in question, incident on the adequate black hole. This can be done in such a 
way because of the fact that Hawking's formula for the emission rate for an outgoing particle at energy oj equals the 
absorption cross section for the same type of particle incoming at energy u. Moreover, outgoing transmission and 
ingoing absorption coefficients are equal. Therefore equilibrium still takes place if the black hole is located in a heat 
bath. 

In what follows we shall investigate massive Dirac fermion emission in the background of (4 + n)-dimensional black 
hole which line element of such a black hole is subject to the relation 

dr"^ 

ds' = ^f{r)de + J^^+ r'dnl^,, (13) 

where / = 1 — I ^o/r j , ro is the radius of the black hole event horizon, while ^ element on 5"+^ 

sphere provided by the relation 

n+2 i~l 

dnl+, = rf^' + E n sin' (14) 

Relation (I12p will constitute the defining equation for massive Dirac fermion fields. By virtue of approximation 
technique we want to achieve the analytical solution of the underlying equation. Namely, we solve the equation for x 
in the near horizon region and then in the far-field limit. Our next task is to match them smoothly. 
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Let us begin with the near horizon limit. After changing of variables in Ea. ((T2)) we arrive at the following: 



dx 
df 



X = 0, 



(15) 



where £ — 2-t2, 

^ n+l 

Then, one redefines xif) = /"(I ~ f)^F{f) and remove singularities at / = and / = 1. The above redefinition 
make it possible to transform this equation to the hypergeometric one of the form as follows: 



dF 



/ (1-/) -J7T + [c-il + a + b)f] —-abF = 



df 



df 



(16) 



It can be checked that the hypergeometric equation parameters satisfy a — a + f3 + £,b = a + P and c — 1 + 2a, while 
a and (3 are given by 



a± 



± 



lujro 
'n+l'' 



(17) 



^± 



where by A we denoted the following: 



(n + 1)2 



-4 



(lu^ - (m - A)2) 
(n+l)2 



(18) 



From this stage on, we shall suppose for simplicity that C{r) = r. Further on, having in mind the criterion for the 
hyperbolic function to be convergent, i.e., Re{c — a — b) > 0, one has to select /3 = /3_. Just the general solution of 
Eq. (|16p may be written in the form 



XNHif) = A_ r (1 -ff Fia + (3 + £,a + P,l + 2a; /) + A+ (1 - FiP + ^ - a, (3 - a,l ~ 2a; /), (19) 



/3 



where A± are arbitrary constants. Because of the fact that no outgoing mode exists near the event horizon of the 
considered black hole, we take a = a- and put equal to zero. This leads to the following solution of equations of 
motion: 



XNHif) = A_ r (1 - F{a + f3 + ^,a + (3,l + 2a; /). 



(20) 



Our next task is to match smoothly the near horizon solution xnh with the far field one in the intermediate zone. 
To do this, first we change the expression of the hypergeometric function near horizon zone from / to (1 — /) by the 
standard relation (see [29| relation 15.3.6). 

Having in mind the limit r — )■ oo or adequately / — > 1 and assuming the low-energy limit u) ^ 1, we are finally left 
with 



XFF = B+ ^/^ JiVT+i>?i^^) + B- Vr Fi yx+4A^('^0: 



(21) 



where YJy and are Bessel functions of second and first order, respectively. 
On expanding in the limit (It — > we conclude that 



XFF{<^r ^ 0) = B_| 



TT 2 



r(i + iVi + 4A2) 

Consequently, taking the near- horizon limit for xnh, one gets 



/I -iVT+iP" , N 



XNHif 



1) - A_ 



-i_iVT+4l^ ^^^i+iVT+Ix^ r(l + 2Q) r(i-e~2/3) 



r(l-e + a-/3) r(l-ha-/3) 



(23) 



^ ^^^yi + iVT+Ix^ ^^^i-i/TW r(l + 2a) r(^ - 1 + 2/3) 



Tii + a + P) Tia + f3) 
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(simplified expression Ea.(|27p) 


I^P (given by Ea.([25|l) 
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0.2 


1.634 X 10"^ 


1.632 X 10"^ 


0.3 


9.287 X 10^^ 


9.244 X 10"^ 


0.4 


3.287 X 10"^ 


3.233 X 10-2 


0.5 


9.004 X 10"^ 


8.612 X 10^2 



TABLE I: Deviation between the values of the absorption probability given by the simplified and complete analytical expression, 
for n = 1, m = 0.01, 5 = 1, Z = and different values of tiro. 

Because of the fact that smooth matching was provided, the ratio of the integration constants may be written in the 
form 



1 r(i + i VI + 4A2) r(i - e - 2/3) r(i VI + 4A2) r(^ + a + /?) r(a + (3) 



B- ^ (ro)^™^ (ici)^™^ r(l-e + a-/3) r(l + «-/3) r(e- 1 + 2/3) 



while the absorption probability is given by the relation 



(24) 



2_ 2»(g* - B) 
~ BB* + i{B*-B) + l ^ ' 

The above relation can be simplified due to the fact that we are using the low-energy limit. In this case BB* ^ 
i{B* — S) 3> 1. Therefore, taking into account the dominant term in the denominator one arrives at the relation 

|.4|'.2.(i-±), (26) 

Using relation ([24|) it can be verified that the absorption probability implies the following: 
I |2_ 8^ ^ (ro)V^+i r(i - (C + /3))2 r(l ~ /3)2 

2^T+4A^ (n + 1) VTTIF r(i - 2/3 - ^2 r( 1 vrT4A2)2 (1 - ^ - 2/3) ■ 

Eq. (j27p is valid for the low-energy range of energy. The low-energy approximation was used during matching the 
two asymptotic solutions in the intermediate zone. However, the simplified analytical relation (j27p is the result of 
series of Gamma function expansions appearing in Eq. ([M|) and therefore its validity is more restricted. In Table 1. 
we presented the values of | ^ |2 derived by using relation ((25|) and equation ((27|) . as lo tq ranges from 0.01 to 0.5, 
for m = 0.01, B = 1, 1 = and n = 1. The results presented in Table 2. are valid for n = 2, the other calculation 
parameters are the same as in Table 1. One can conclude that, the agreement between these two values in question 
is remarkable as Q ro reaches the value 0.5, then the deviation between these values appears. 

We remark that expression (|27p is quite general, valid for both bulk and brane massive Dirac fields. The case of 
bulk fermions is connected with the eigenvalues of Dirac operator on a transverse manifold which is S^''^'^ sphere. In 
case under investigation one should consider (n + 2)-dimensional sphere. The eigenvalues for spinor where found 
in Ref. (s^l • They imply the following relation: 



where I = 0, 1, . . . On the other hand, massive brane fermions, live on four-dimensional brane and propagate in the 
gravitational background provided by 

r7r2 

ds^ = -f{r)de + J^^+ r^dn", (29) 
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Loro \A\'^ (simplified expression Eq.^^) \A\'^ (given by Ea. pS]) ) 

0.01 2.153 X 10"" 2.153 x 10"" 

0.05 5.921 X 10"* 5.921 x 10"* 

0.1 2.041 X 10"^ 2.041 X lO"** 

0.2 7.186 X 10"^ 7.186 x lO""^ 

0.3 5.865 X 10"* 5.863 x 10"'' 

0.4 2.641 X 10"^ 2.638 x 10"^ 

0.5 8.608 x 10"^ 8.571 x 10"^ 



TABLE II: Deviation between the values of the absorption probability given by the simplified and complete analytical expression, 
for n = 2. Other calculation parameters are the same as in Table 1. 



n+l 



where / = 1 — ( f o/r j , ?'o is the radius of the black hole event horizon. This metric on the brane [J [3| arises 



when we project out all angular variables that parameterize the extra dimensions, i.e., if we set tpi = 7r/2 for i > 2. In 
the above line element n stands for the number of extra dimensions that can exist transverse to the considered brane. 
Then, the eigenvalues on a transverse manifold are of the form — (l + 1)^. 

The other kind of n-dimensional black hole which can be considered by relation ([27|) is a tense brane black hole. 
It turned out that the most examinations of extra-dimensions black hole and their physics were devoted to the zero 
brane tension case. In principle finite brane tension ought to modify the physics of a kind of black hole. The nonzero 
tension on the brane can curve the brane as well as the bulk. It was shown that a tense brane black hole is locally a 
higher-dimensional Schwarzschild solution 3l| pierced by a tensional brane. This caused that a deficit angle appeared 
in the {n + 2)-dimensional unit sphere line element. 

Some attempts to examine the problem in question were conducted. Namely, studies of massless fermion excitation 
on a tensional three- brane were carried in Ref.jsll, while the late-time behaviour of massive scalar hair in the back- 
ground of n-dimensional tense brane black hole was studied in (33j . Examinations of Hawking radiation of massive 
scalar fields in the background of a tense five and six-dimensional black hole were carried in Ref. 34]. On the other 
hand, emissions of massless scalar fields into the bulk from six-dimensional rotating black hole pierced by a three- 
brane were studied in [s^. Ref. [sB] was devoted to the numerical studies of evaporation of massless scalar, vector and 
graviton fields in the background of a six-dimensional tense brane black hole. Growing interests in codimensional-2 
braneworlds lead us to modify the gravitational action by implementing Gauss-Bonnet term or to consider black hole 
solutions on a thin three-brane of codimension-2 [s^, HI] . 

In the case of a tense brane black hole 5"+^ sphere is threaded by a codimension-2 brane, so the range of the one 
of the angles, let us say, ipi will he < < 2ttB. Parameter B measures the deficit angle about axis parallel with 
the brane intersecting the sphere in question. For such a kind of transverse manifold the eigenvalue for spinors will 
be of the form uM 



A = /c H h c 



(30) 



where c = ±1/2, ±3/2, .... 

To complete this section let us turn to study the luminosity of the Hawking radiation for massive Dirac field. It is 
provided by the expression 



L = 



e^BH -j- 1 



(31) 



4(ro)™^+i r(l - (C + /3))2 r(l - /3)2 



2^/l+4A^ (n + 1) VI + 4A2 r(l - 2/3 - 0^ r(iVl + 4A2)2 (1 - ^ - 2/3) 
1 (TBff)^/^+^+' C(Vl + 4A2 + 3) r(v/l + 4A2 + 3) 



02 m 



V1 + 4A2 (rsH)^^+^+' C(Vl + 4A2 + l) r(Vl + 4A2 + l) 



7 



where Tbh — '^j^^ is the temperature of the black hole while C, is Riemann zeta function, ai and 02 are given by 

oi = 1 ]- , a2 = 1 7^ ■ (32) 

In the above relation pip by A we mean A |/=o- In our low-energy approximation the luminosity of the black hole 
Hawking radiation is calculated with respect to ? = mode. However, for the sake of completeness we write in Eg. ([51]) 
the integral range from zero to infinity. One should have in mind that our analysis has focused on the low-energy 
spectrum and the value of luminosity is based on the lower part of the spectrum. It can happen that modifications 
may appear for high-energy part. 



IV. THE ABSORPTION PROBABILITY AND HAWKING RADIATION IN THE SPACETIME OF 

HIGHER DIMENSIONAL BLACK HOLE 

In this section we present plots of the absorption probability and luminosity of Hawking radiation for different 
parameters in different kinds of n-dimensional black holes. The absorption probability is a dimensionless constant 
and it should range from to 1 for the whole energy range. Our considerations are confined to the low-energy limit 
therefore our plots will cover only a part of the whole plot for | A p. The whole plot for | .4 ^ can be found by 
numerical calculations (see, e.g., Fig.l in Ref.flB] for the comparison of absorption probability obtained in the low- 
energy limit and those get by various numerical approximation schemes). 

In all our logarithmic plots we use base- 10 units. One also elaborates the dependence of the absorption probability 
and luminosity of Hawking radiation for massive hulk and hrane fermions. 



A. Bulk Emission 



We begin our considerations with studies of massive hulk Dirac fermions emission. In Fig.l we plotted \ A]"^ for 
massive fermion fields with respect to lo for different number of extra dimensions n = 1 , ... 4 in the background 
of (n -f 4)-dimensional Schwarzschild black holes and in the spacetime of tensional brane black hole. We fixed the 
multipole number I = for the Schwarzschild case, while for the tensional brane we put fc = 0, c — 1/2 and B — 0.8. 
Other calculation parameters are: m = 0.01 and tq = 1. One should recall that the event horizon radius of a tense 
brane black hole is connected with the Schwarzschild radius tq by the following relation: 

rTBBH = — (33) 
B— 

It turned out that the absorption probability for massive hulk Dirac fermions for the considered tense brane black 
hole is smaller comparing to the (n + 4)-dimensional Schwarzshild case. 

It was revealed in Ref. _16J using the WKBJ approximation and the Unruh method that the emission rate of the bulk 
massless fermions was dependent on the the dimensionality of the spacetime. The bigger is the spacetime dimension 
the smaller emission rate we get. It turned out that at a certain intermediate energy there was a region where the 
emission rate became approximately independent of the dimension of underlying background. Our Fig.l confirms this 
tendency, i.e., for massive Dirac fermions the absorption probability decreases with the increase of dimension of the 
spacetime. 

In Fig. 2 we present the absorption probability versus lo for n — 5, and 6, for different multipole numbers: / = 0, 1, 
and 2, in the Schwarzschild background. We also examine the case of a tense brane black hole and put fc = 0, 1, 
and 2. Other calculation parameters we take into account are: m = 0.01, B = 0.8, rg = 1, and c = 1/2. One can 
notice that the absorption probability decreases as the multipole number increases. 

In Fig. 3 we depict | ^ ^ versus uj for different mass of the Dirac massive fermion fields in five and six-dimensional 
Schwarzschild black holes background. We consider the case when m = 0,01, 0.1, and 0.15. At the beginning | A P 
for massive Dirac fermions decreased as the mass of the field increased (the same situation was also revealed in studies 
of the bulk absorption probability for scalars in the spacetime of Schwarzschild black hole located on a three-brane of 
finite tension (see Fig. 4 in Ref.[l3])) but then there was a region where the inverse situation took place. Contrary to 
the previuos behaviour the absorption probability increased with the growth of the particles masses. This behaviour 
may be seen on the left panel for five-dimensional spacetime and the same tendency occurs in six-dimensional case. 

The dependence of the absorption probability on the Schwarzschild black hole event horizon radius is presented in 
Fig. 4. We studied the case of tq = 1.2, 1, and 0.8. One can find that the absorption probability increases with the 
growth of radius of the black hole in question. 
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On the other hand, inspection of Eq. (|27|) reveals that the absorption probabiHty for massive bulk fermions is strictly 
bounded with the parameter B characterizing the brane tension. In Fig. 5 we plotted | A P versus uj for different 
spacetime dimensionality n = 5, and 6 and for different values of parameter B. In the considerations we put 
B = 1, 0.9, 0.8, respectively. Fig. 5 provides the fact that the smaller B (brane tension increases) the smaller 
absorption probability one gets. 

In Ref.[35| the emission bulk massless scalar fields from a six-dimensional rotating black hole pierced by a three- 
brane was numerically studied. For the low-energy region the authors find the dependence of absorption probability 
on parameter B. The curves presented in Fig. 2 in (35j for the case of angular parameter equal to zero has the same 
tendency as our Fig. 5 (right panel). The same situation takes place when one considers massive scalar fields emitted 
by a tense brane black hole [13] ■ So one can conclude that the brane tension suppresses the emission of massless 
scalars, massive scalars as well as massive bulk fermions. 

In Fig. 6 we present the | A P for five and six-dimensional tensional brane black holes for different c = 1/2, 3/2, 5/2, 
respectively. Other calculation parameters are: m = 0.01, ro = 1, A; = and B ~ 0.9. It occurs that the bigger c is 
the smaller absorption probability one obtains. 

B. Emission on the Brane 

Now, we proceed to study the absorption probability of Dirac massive fermions propagating on the brane. In the 
left panel of Fig. 7 we present the brane absorption probability for fixed I = for different spacetime dimensionalities 
n = 1,...4. It turns out that the absorption probability increases with the increase of spacetime dimensionality 
n for the brane black hole which all angular variables parameterizing the extra dimensions are projected out. The 
same behaviour of | ^ ^ was shown for massless brane fermions ^ in the spacetime of rotating brane black hole. 
These results were also confirmed numerically in Ref.[l3| where it was also revealed that any increase of spacetime 
dimensionality or angular parameter of black hole enhances the massless fermion emission rate of the considered black 
hole. Studies of radial Teukolsky Eq. conducted in Ref.jl^ also confirmed the tendency of behaviour we obtained 
(see Fig.l from Ref.fisj for the rotation parameter equal zero depicting greybody factor for brane localized fermions). 

In the right panel of Fig. 7 we plotted the | p as a function of oj for different multipole number Z = 0, 1, and 2 
in five and six-dimensional spacetimes. It happens that the brane absorption probability decreases as the multipole 
number increases, as was the case for the bulk radiation. 

Next, in Fig. 8 we examined the behaviour of the absorption probability on the brane for different masses of the 
fermion particles. We put m = 0.01, 0.1, 0.15, respectively. At the beginning | A P decreases as the mass of fermions 
increases but there is a region where the situation changes, i.e., for the increase of mass one gets also increase of the 
absorption probability. The same situation was observed in the bulk case. 

On the other hand, in Fig. 9 we plotted the dependence of | ^ p on the radii of black hole event horizons. One 
studies the cases of ro = 1.2, 1, and 0.8. It can be concluded that the absorption probability for massive brane Dirac 
fermions increases as the event horizon radius increases. The same behaviour we have in the case of massive bulk 
fermions. 



C. Luminosity 

Finally, we plotted the luminosity of the Hawking radiation L as a function of m for different number of extra 
spacetime dimensions n = 1, ... 4. We plot the luminosity of the Hawking radiation for the mode I — which plays 
the dominant role in a greybody factor in the low-energy approximation. In the left panel of Fig. 10 we present 
the bulk absorption probability for Schwarzschild and tensional brane black holes. In the right one, we plotted the 
brane absorption probability in the background of higher dimensional Schwarzschild black holes. One can remark 
that the luminosity of Hawking radiation increases with the increase of n, spacetime dimensionality for massive bulk 
and brane fermion emission. The other salient feature is that luminosity for brane massive fermion is substantially 
higher comparing to the luminosity for bulk fermions for the corresponding dimensionality of spacetime. Of course 
this behaviour is valid only for the low-energy limit. The behaviour of the luminosity valid for the whole energy range 
should be obtained by numerical studies. However our analytical results tell us about the tendency of behaviour of the 
luminosity in this energy limit. In Ref. [lB| the massless fermion emission rate was studied numerically and it was shown 
that for the brane localized emission one obtained the same tendency of behaviour as ours, i.e.. Hawking radiation 
increased with the growth of the dimensionality of the spacetime. It was also found that in the intermediate energy 
range the emission rate became independent on the dimensionality. However this case is beyond our approximation 
scheme. 
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The same tendency of behaviour was revealed in Ref. Ts*] for the massless scalar emission rate in the spacetime of 
n-dimensional rotating black hole. Namely, the authors find that for the rotational parameter put to zero the total 
emission depends strongly on the number of dimensions both for bulk and brane emission. It turned out that on 
average, more energetic emission is in the bulk than on the brane. Unfortunately, our plots of the luminosity for 
massive Dirac fermions can not confirm this behaviour valid for massless scalars because of the considered energy 
limit. 

The dependence of a flux emission spectra of massless fermion on the brane emitted from rotating black hole on the 
spacetime dimension was also confirmed in Ref. . 

Thus, summing it all up we can conclude that in the considered energy limit our key result is that the luminosity 
for massive Dirac fermions strongly depends on the bulk and brane dimensions. But as one can try to consider the 
whole energy range these results should be taken with a large grain of salt. 



V. CONCLUSIONS 



In this paper we studied Hawking emission of massive Dirac fermion fields in the spacetime of static (4 + n)- 
dimensional black holes. We elaborated the case of (4 + n) Schwarzschild black hole, (4 + n) tense brane black hole as 
well as propagation of brane fermions in the spacetime of brane black hole when all angular variables parameterizing 
the extra dimensions in the considered line element were projected out. It happened that the treatment of Dirac 
fermions in spherically symmetric spacetime was simplified to great extent due to the few properties of the Dirac 
operator. The crucial quantities characterizing emission of massive Dirac fields are dependent on the eigenvalues A 
of the Dirac operator on the so-called transverse manifolds. We have elaborated analytically equations of motion 
for the aforementioned degrees of freedom and by means of matching technique one finds in the low-energy limit 
an analytical expression for the absorption cross section, luminosity of Hawking radiation. We derive quite general 
formulae for these quantities which enable one to treat all these spacetimes in question. Our analytical considerations 
are supplemented by plots of | p and the luminosity of Hawking radiation L. We found that | ^ p for bulk massive 
Dirac fermions in the spacetime of (n -|- 4)-dimensional Schwarzschild black hole and in the background of a tense 
brane black hole decreases with increasing of the dimensionality of the considered spacetime. In higher dimensional 
Schwarzschild spacetime the absorption probability decreases as the multiple number I increases. For this spacetime 
it was also revealed that | A P depended on mass of the emitted field. First, at the beginning of the energy range, we 
have the situation that the bigger mass of Dirac fermion is the smaller | -4 p one gets. Then, this tendency changes 
and we observe the increase of the absorption probability with the growth of mass of particle in question. It was 
also observed that the absorption probability increased with the increase of the radius of the event horizon. We note 
that I ^ P reveals its dependence on parameter characterizing tense brane black hole. Namely, the smaller B one 
considers the the smaller \ A]"^ one gets. One ought to have in mind that B is bound with the brane tension, i.e., the 
smaller B the greater tension is exerted on the brane black hole. In the case of a tense brane black hole the growth 
of a multipole number c implies the decreasing of the absorption probability. 

We also studied the case of brane massive Dirac fermion field. It was shown that | ^1 P increased with the increase 
of n for brane black hole which all angular variables parameterizing the extra dimensions were projected out. The 
same character of behaviour was shown in Ref. 'i'l for the case of massless brane fermions. It was also noticed that 
\ A\'^ for different masses of brane fermions behaved in the similar way as in the case of bulk massive Dirac fermions. 
On the other hand, | A P analyzed for different event horizon radius also reacts in the similar way as in the case 
of bulk fermions. Finally, analyzing luminosity of Hawking radiation for massive brane Dirac fermions it was found 
that the brane luminosity was substantially higher comparing to the luminosity of Hawking radiation for bulk massive 
fermions for the corresponding values of spacetime dimensionality n. One should remark that the behaviour of the 
luminosity is valid for the low-energy limit, when li tq <C 1. 
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FIG. 1: Absorption probability | ^ |^ for bulk Dirac fermions with mass m as a function of w, for different space dimensionality 
n = 1, 2, 3, 4 (curves from the top to the bottom, respectively). The solid lines represent the case of (n + 4) -dimensional 
Schwarzschild black holes, the dashed lines represent (n+4)-dimensional tense brane black holes with B = 0.8. Other calculation 
parameters are: c = 1/2, fe = 0, 1 = 0, ro = 1, and m = 0.01. 




FIG. 2: Absorption probability | A |^ for massive Dirac field as a function of co, for difi'erent values of / = 0, 1, 2 (solid 
lines, curves from the top to the bottom, respectively) for Schwarzschild black holes and for fc = 0, 1, 2 for tense brane black 
holes (dashed lines, curves from the top to the bottom, respectively), for five and six-dimensional cases (left and right panel, 
respectively). The rest of the calculation parameters are the same as in Fig.l. 
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FIG. 3; Absorption probability | A for massive formion particles in a spacetimo of five and six-dimensional Schwarzshild 
black hole (left and right panel, respectively), for different values of the mass of Dirac fermions m = 0.01, 0.1, 0.15 (curves 
from the top to the bottom, respectively). Other calculation parameters are: I = 0, and ro = 1. 




FIG. 4: Absorption probability | .4 |^ for massive fermion particles in a spacetime of five and six-dimensional Schwarzshild 
black hole (left and right panel, respectively). Curves from the top to the bottom are for different values of black hole event 
horizon radius ro = 1.2, 1, and 0.8. Calculation parameters are: I = 0, and m = 0.01. 




FIG. 5: Absorption probability | A f versus ui for different values of -B = 1, 0.9, 0.8 (curves from the top to the bottom, 
respectively) and for different spacetime dimensionality n = 1, and 2 (panels from left to the right) for tense brane black holes. 
The rest of the parameters are the same as in Fig.l 
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FIG. 6: Absorption probability | A P versus lo for massive fermiori field for different values of c = 1/2, 3/2, and 5/2 (curves 
from the top to the bottom, respectively) and for different spacetime dimensionalities n = 1, 2 (panels from left to the right) 
in the background of a tense brane black hole. Other calculation parameters are: m = 0.01, ro = 1, k = 0, and B = 0.9. 




FIG. 7: Left panel: absorption probability | A \'^ versus lo for massive branc fcrrnions for fixed / — and different spacetime 
dimensionalities n = 1, 2, 3, 4 (curves from the bottom to the top, respectively). Right panel: absorption probability for 
different multipole numbers 1 = 0, 1, 2 (curves from the top to the bottom, respectively) for five and six-dimensional cases 
(dashed and solid lines, respectively). Other calculation parameters: m = 0.01, and ro = 1. 




FIG. 8: Absorption probability | A |^ for massive brane fermion particles in a spacetime of five and six-dimensional Schwarzschild 
black holes (left and right panel, respectively), for different values of m = 0.01, 0.1, 0.15 (curves from the top to the bottom, 
respectively). Other calculation parameters are: I = 0, and ro = 1. 
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FIG. 9: Absorption probability | A P for massive brane fermion particles in a spacetime of five and six-dimensional Scliwarzschild 
black holes (left and right panel, respectively). Curves from the top to the bottom are for different values of black hole event 
horizon radius ro = 1.2, 1, and 0.8. Other calculation parameters are: I = 0, and m = 0.01. 
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FIG. 10: The luminosity of the Hawking fermion radiation L versus rn propagating in the bulk (left panel) on a Scliwarzschild 
black hole spacetime (solid lines) and on a tense brane black hole spacetime with B = 0.9 and fc = 0, c — 1/2 (dashed lines) for 
different space dimensionalities n = 1, 2, 3, 4 (curves from the bottom to the top, respectively). Right panel: the luminosity 
of Hawking radiation L versus m propagating on the brane for different space dimensionalities n = 1, 2, 3, 4 (curves from the 
bottom to the top, respectively). For all cases we put ro = 1. 



